Slowing and stopping of chemical waves in a narrowing canal 
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The propagation of a chemical wave in a narrow, cone-shaped glass capillary was investigated. 
When a chemical wave propagates from the wider end to the narrower end, it slows, stops, and then 
disappears. A phenomenological model that considers the surface effect of the glass is proposed, 
and this model reproduces the experimental trends. 
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INTRODUCTION 

The mechanism of information processing in neuronal 
systems is one of the central issues in modern science, 
including biology, biochemistry, and physics. It is known 
that information processing in nerve cells involves a tem- 
poral change in the membrane potential. It has been 
well-established that, in cablelike nerve cells, pulses on 
the membrane propagate at a constant velocity and am- 
plitude P, 0] . Recently, on the basis of advances in ex- 
perimental techniques, it is becoming clearer that the 
manner of this pulse propagation is critically dependent 
on the thickness or width of the nerve cells. For exam- 
ple, it has been reported that the dendritic shape of nerve 
cells strongly affects the propagation of an excited wave 

^ IB| ' These effects enable information processing em- 
bedded in time-dependent signals, such as coincidence 
detection, in which the timing of a pair of inputs can be 
determined [lllli. 

On the other hand, it has been claimed that the essence 
of the neuronal system can be understood in terms of a 
simple differential equation, the FitzHugh-Nagumo equa- 
tion, as a type of reaction-diffusion system |2lUi|. Thus, 
a neuronal pulse is a manifestation of a dissipative struc- 
ture generated in nonequilibrium open systems . Sev- 
eral experimental and theoretical studies have been un- 
dertaken on the reaction-diffusion systems as a model of 
the neuronal system. 

In such studies, the Belousov-Zhabotinsky (BZ) re- 
action is often used as an experimental model of the 
reaction-diffusion systems In this system, chemi- 

cal waves propagate at a constant velocity and amplitude 
and, thus, exhibit the characteristics similar to those of 
pulses in nerve cells. The characteristics of the BZ re- 
action can be well described using a numerical model, 
the Oregonator [13. For example, the Oregonator 

model reproduces the existence of oscillatory and ex- 
citable states, the formation of fascinating patterns (such 



as target patterns and spiral patterns) fl^ , and also fea- 
tures associated with hydrodynamical effects The 
Oregonator model has a mathematical structure that is 
similar to the FitzHugh-Nagumo equation. Therefore, 
to better understand information processing in the neu- 
ronal system, it may be useful to compare the behav- 
ior of chemical waves in the BZ reaction, although there 
exists a fundamental difference between them; a neu- 
ronal pulse propagates on the two-dimensional tubular 
membrane and a chemical wave propagates in the three- 
dimensional volume. 

Over the past decade, the behavior of the chemical 
waves in a narrow space has been studied using a gel 
system [tH. 0, droplet system J^], and bead system 
[20I |2]| . Toth and co-workers reported that a chemi- 
cal wave does not propagate when it expands from a thin 
capillary to the bulk solution. They claimed that a chem- 
ical wave can perform information processing |22l l23l |. 
Masere et al. performed experiments on the propagation 
of a chemical wave in a vertical capillary. They discussed 
their experimental results with special emphasis on the 
gravitational effect ,2^]. In these previous studies, the 
features of the chemical-wave propagation in a narrow 
straight capillary with a constant diameter were exam- 
ined. 

In the present study, we show that a chemical wave 
slows, stops, and then disappears in a narrowing glass 
capillary. We discuss this phenomenon in term of the 
surface volume ratio for the reaction field, i.e., the surface 
can seriously affect propagation of the chemical wave. 



EXPERIMENTS 

All chemicals were analytical-grade reagents and used 
without further purification. An aqueous solution of 
ferroin, or tris (1,10-phenanthroline) iron (II) sulfate, 
was prepared by mixing stoichiometric amounts of 1,10- 
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FIG. 1: Schematic representation of the experimental system; 
the capillary was arranged horizontally and observed from 
above by a charge-coupled device (CCD) camera. 

phenanthroline and ferrous sulfate in pure water. The 
water was purified with a Millipore MiUi-Q system. 
The BZ medium contained 0.35 M sodium bromate 
(NaBrOa), 0.39 M sulfuric acid (H2SO4), 0.12 M malonic 
acid (CH2(COOH)2), 0.05 M sodium bromide (NaBr), 
and 4.0 mM ferroin (Fe(phen)3''") (excitable condition). 
We allowed the BZ solution to stand for ~ 10 min and 
then stirred it to diminish the formation of bubbles in 
the solution. 

We prepared a glass capillary (borosilicate glass with 
a length of 100 mm and inner radius of 0.75 mm; Sut- 
ter Instrument Co.) with a decreasing width by pulling 
a straight capillary with a Micropipete puller (Model P- 
97/IVF, Sutter Instrument Co.). The both ends of the 
capillary were open, and the gradient of the inner radius 
was 0.036. To make it easier to show the experimental 
results, the x-axis is set as shown in Fig. ^ The capil- 
lary was immersed in the BZ medium, fulfilled with the 
BZ medium inside it, and then situated on a horizontal 
plate for observation with a charge- coupled device (CCD) 
camera (Digital Microscope Unit, Keyence). As soon as 
bubbles formed inside the capillary, we stopped the ob- 
servation: the data reported in the present study are only 
from the experimental series without bubble formation. 
All the experiments were performed at room tempera- 
ture (20 ±3 ° C). Images were recorded on a videotape 
and analyzed by an image processing system (Himawari, 
Library, Inc.). 



RESULTS 

A chemical wave was initiated at the wider end by a 
silver wire without disturbance. A chemical wave prop- 
agating toward the narrower end slows, stops, and even- 
tually disappears. The oxidized solution in the capillary 
then reverts to the reduced state. A spatio-temporal plot 
of images of the capillary along the long axis is shown in 
Fig. |2Ia), and the snapshots (taken every 5 s) are shown 
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FIG. 2: Experimental results on the propagation failure in 
a narrow capillary (the chemical wave propagated from left 
(wider end) to right (narrower end)): (a) Spatio-temporal 
plot of chemical-wave propagation and (b) snapshots with an 
interval of 5 s. When the wave propagated to the narrower 
region, it slowed, stopped, and then disappeared. Darker and 
brighter areas correspond to reduced and oxidized states, re- 
spectively. 

in Fig. |2Ib). These clearly show that the velocity of 
the chemical wave gradually decreases before it stops. 
This feature was observed in every chemical wave initi- 
ated. With each subsequent wave, the position at which 
the chemical wave stops has a tendency to shift slightly 
toward the narrower end. Figure |31 shows the tempo- 
ral change in the position (Fig. Ela)) and velocity (Fig. 
Efb)) of the chemical wave, and the velocity depending 
on the position (Fig. EIc)). From the observation, the 
point where the chemical wave propagation stopped Xg 
was 2.7 ± 0.2 mm. We have confirmed the reproducibil- 
ity of the experimental trend through a number of several 
repeated experiments. 

DISCUSSION 

Here, we discuss here the characteristic change in the 
manner of wave propagation in a glass capillary. Accord- 
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FIG. 3: Features of the propagation failure of chemical waves 
in the experiments as shown in Fig. |21 (a)temporal change 
in the position of the chemical waves, (b)temporal change in 
the velocity of the chemical waves, and (c)the relationship 
between position and velocity. 

ing to an analysis with the Oregonator, the velocity of 
the chemical wave in the BZ reaction can be written as 

m 

V « i ihHADf' , (1) 

where H is the concentration of H+ ions, A the con- 
centration of bromate ions, D the diffusion constant of 
the activator, and the rate constant of the compo- 
nent reaction, which is introduced as an autocatalytic 
process by summation of the two elemental processes in 
the Oregonator model [2^. Equation |^ is valid for a 
plane and steady traveling wave but does not apply to 
a wave that shows a change in velocity. By assuming 
that the concentration dependence of the local velocity 
along the capillary is also given in a similar way, we sim- 
ply adapt eq.Q) to interpret the experimental trends in 
a semiquantitative manner. 

The X-axis is set in the direction of wave propagation. 
We describe the inner radius of the glass capillary (r) by 
assuming it to be a part of a cone: 

r = a{xo - x), (2) 
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FIG. 4: Plot of vs 1 — b/[a(xo — x)]. Five chemical waves 
were measured, and the different marks correspond to the 
different waves. By line-fitting, the gradient was calculated 
to be 450 mm"^ s^. 
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FIG. 5: Relationship between the position and velocity of 
a chemical wave, based on the theoretical discussion given 
in eq.ljlj; the parameters used here are determined from the 
experiments. 



where a corresponds to the gradient of the inner radius 
of the glass capillary and xq is the position that would 
be the top of the cone. 

The concentration of hydrogen ions (H) can be de- 
creased by the effect of the glass surface of the capillary. 
Note that H is the average concentration over the cross 
section, and that the effect on the hydrogen concentra- 
tion by the glass surface is effective only near the surface. 
Amemiya et al. reported that the effective concentration 
of hydrogen ions is not homogeneous in the porous glass 
system |2a|- Thus, the Na'^ ions that are connected to 
the Si02~ groups could be replaced by H+ ions, which 
causes a decrease in the value of H near the glass surface. 
We assume that this effect is proportional to the surface 
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volume ratio: 



2'Kr 



H I 



(3) 



where 6 is a positive constant that reflects the effect of 
the glass surface. Although, at present, it seems difhcult 
to evaluate b experimentally, as a future subject, it may 
be valuable to try to measure the change in the chemical 
environment near the surface. The parameter b (= 26) 
is a positive constant that is introduced to simplify the 
equation. 

From the abovementioned equations, we can derive the 
relationship between the position and velocity: 



V = Vq 



a{xo — x) 



1/2 



(4) 



where vq is the velocity of the chemical wave without 
any effects from the glass surface. By letting = in 
eq.Q, we can determine the position Xg, which is where 
the chemical wave stops: 



h 

Xs= Xq . 

a 

Equation Q can also be written as 

Vq) a{xo-x)' 



(5) 



(6) 



From the experimental results, the parameters were de- 
termined to be Xs = 2.7 ± 0.2 mm, Xq = 5.5 mm, and 
a = 0.036. Using these values, we plotted the relation- 
ship between v"^ and 1 — b/{a{xo — x)} in Fig. 0] where 
b is calculated as b = 0.20 ± 0.01 using eq. ©. The 
results show a linear relationship between these two val- 
ues. By line-fitting, the gradient was calculated to be 
450 ± 50 mm~^ s^. This value corresponds to I/vq'^, in- 
dicating that Vq — 0.047 ± 0.003 mm/s, which roughly 
corresponds to the wave velocity obtained through the 
experiment using a wide capillary (0.06 ± 0.01 mm/s). 
Using this value, the plot between v and x is shown in 
Fig. |S1 The relationship between v and x in the exper- 
iments, as shown in Fig. |2tc), was well-reproduced in a 
quantitative manner. The velocity profile near the point 
where a chemical wave stops can be thus explained. 

In the abovementioned discussion, the decrease in the 
effective concentration of the hydrogen ion, (i?cff), near 
the glass surface seems to be induced by the adsorption 
by the Si02~ groups on the glass surface. With time, this 
substitution continues, and the effect of the glass surface 
on H decreases. This may be the cause of the slight shift 
in the position at which the chemical wave stops. 

Because we performed the experiments under an oxy- 
gen atmosphere, it is possible that the oxygen disturbs 
the propagation of the chemical wave. In regard to this 



problem, we conducted the observation using a straight 
capillary with open ends and confirmed that absorption 
of the oxygen from an open end has a negligible effect on 
the manner of wave propagation. Therefore, it is most 
probable to expect that the significant factor is the ab- 
sorption of the H+ ion to the inner wall of the glass cap- 
illary. 



CONCLUSION 

When a chemical wave propagates from the wider end 
of a glass capillary to the narrower end, it slows, stops, 
and then disappears. This failure to propagate can be 
caused by the effect of the glass surface inside the cap- 
illary. A phenomenological model was considered, and 
the essential aspects of the slowing and stopping of the 
chemical wave were reproduced. 

This study may promote a further understanding of 
signal transduction in the neuronal system, especially in 
narrow tubular parts, such as in the dendrites of nerve 
cells. 
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